ABSOLUTE RUIN IN THE ORNSTEIN-UHLENBECK TYPE RISK MODEL 
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Abstract. We start by showing that the finite-time absolute ruin probability in the classical 
risk model with constant interest force can be expressed in terms of the transition probability 
of a positive Ornstein-Uhlenbeck type process, say X. Our methodology applies to the case 
when the dynamics of the aggregate claims process is a subordinator. From this expression, we 
easily deduce necessary and sufficient conditions for the infinite-time absolute ruin to occur. We 
proceed by showing that, under some technical conditions, the transition density of X admits 
a spectral type representation involving merely the limiting distribution of the process. As a 
by product, we obtain a series expansions for the finite-time absolute ruin probability. On the 
way, we also derive, for the aforementioned risk process, the Laplace transform of the first-exit 
time from an interval from above. Finally, we illustrate our results by detailing some examples. 
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1. Introduction 

Let Z = > 0) be a driftless subordinator defined on a filtered probability space 

(0, (J-t)f>o, P), that is Z is a M"*"— valued process with stationary and independent incre- 
ments. It is well known that the law of Z is characterized by its Laplace exponent which admits 
the following Levy-Khintchine representation 

(1.1) </.(/3) = -logE(e-''^^)= / {l-e-Py)u{dy), /3 > 0, 







where the Levy measure satisfies the integrability condition /q°°(1 A y)u{dy) < oo and we refer 
to the monograph of Kyprianou |20] for background on subordinators. By P we denote the law 



of the dual of Z, i.e. the law of —Z under P. For x,r,c G M, we denote by pi^''^^ the law of the 
process X defined, for any t > 0, by 

(L2) Xt = e^* (^x + e-^M(cs - Z,) 

where Z ~ P. Similarly, Fx'"^"^ stands for the law of the process X as defined in ()1.2p with Z ~ P. 

We simply write (resp. Pi^^) for P^'"''^^ (resp. pi'"''^^). Note that equivalently, X is the unique 
strong solution to the stochastic differential equation 

(L3) dXt = {rXt + c)dt - dZs, Xo = x. 
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The process X under ' with c > has been used in the hterature to model the reserves 
of an insurance company, the parameter c standing for the premium rate, the jumps of Z 
standing for the claims and x standing for the initial value of the reserves. In particular, when 
z^(0, c«) < oo, X under P^''''^^ is the classical risk process where the claims arrive according 
to a Poisson process with intensity parameter z^(0, oo) and the claim distribution is given by 

z^(dx)/i^(0, oo). When r > 0, the process X under P^'^''^^ has been suggested as a risk process 
where the cost of lending/borrowing money are taken into account. In this model the insurer 
earns (credit) interest when the surplus is positive and when the surplus becomes negative, the 
insurer can cover the deficit by a loan for which he has to pay a (debit) interest. Although in 
practice the debit interest is much higher than the credit interest, we restrict ourselves to the 
case where both rates are equal since this choice is particularly tractable and allows us to use 
techniques which can no longer be used in the general case. Processes of the form (II. ID with Z 
a general Levy process are known in the literature as processes of Ornstein-Uhlenbeck type (for 
short OU-type) and therefore in this paper we call the process X under P^^''^^ with r > 0, the 
OU-type risk process. 

From (jl.3p . it is clear that when the OU-type risk process reaches the interval (—oo,— c/r], 
the premium rate c can no longer compensate the interest payments and so the surplus will 
decrease to minus infinity. Following Gerber [15], we say that in this case absolute ruin occurs. 
This model and in particular the event of absolute ruin, has been the focus of many research in 
insurance mathematics, its first appearance can be traced back to Segerdahl [35j. For more recent 
investigations and substantial refinements, we mention Gerber [15], Dassios and Embrechts [10] . 
Embrechts and Schmidli [T2], Schmidli [33], Sundt and Teugels [37], Albrecher et al. |2], Gerber 
and Yang [IT] and Cai We refer to the survey papers of Paulsen [28] and [27] for an overview 
of ruin models with interest. For general background in ruin theory, we refer to Gerber [16] and 
Albrecher and Asmussen [1]. 

We also point out that the process X under Wx with r > 0, is well known in the literature 
and has appeared in various settings. For instance, in mathematical finance, this process has 
been used by Barndorff-Nielsen and Shephard for the modeling of the stochastic volatility of 
stock prices, see j3j. It also belongs to the class of one factor affine term structure models, 
see e.g. Filipovic [T3|. Moreover, when the Levy measure is finite, it is a specific example of a 
Poisson shot noise process, see e.g. Perry et al. |29j and Iksanov and Jurek [19j. 

The remainder of the paper is organized as follows. The next section is devoted to the 
statement of our main results which are then proved in Section HI Examples illustrating our 
approach are presented in Section [3l 

2. Main results 

We start by providing a representation of the law of the absolute ruin time. To this end, we 
introduce some notation. First, let, for any a < x, 

Ta = inf{s > : Xg < a] 

be the first-passage time below the level a for X. Henceforth, we shall assume that r > 0. We 
are interested in the distribution of t_^/j, under P^'^''^-*. We call this random variable the absolute 
ruin time. As mentioned in the introduction, the reason for the adjective 'absolute' is that once 
the OU-type risk process goes below the level —c/r, it will never go back above this critical 
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level, i.e. 

(2.1) p(?-'^)(Xt>-c/r,t>T_,/,)=0. 

We first note that from (|1.2p we immediately see that the process X under P^*^'^^ has the same 
law as the process {Xt — c/r,t > 0) under P^'^^y^. In particular 

Pi^''^) (t_,/, G dt) = P^,/, (to G dt) , t>0. 

Based on this observation, we state all the results in the paper for the c = case only; the 
analogue for c 7^ is then obvious. The first theorem gives the link between the distribution of 

(r) r) 

To under P^ ' and the distribution of X under Pq ^ , which leads to an explicit expression for 
the Laplace transform in space of the finite-time absolute ruin probability. This is in contrast 
with the the finite-time ruin probability in the r = case, where only an explicit expression for 
the double Laplace transform in space and time exists (cf. Theorem 8.1(ii) of Kyprianou [20J). 

Theorem 2.1. For any a; > and t > 0, we have 

(to <t) = P-'^ (Xt > x) . 
Consequently, for any /3 > 0, we have 

P5(r0 >t) = exp 1^ 0(/3e-^'^)ds^ 

where stands for the exponential distribution of parameter /3 > (where we used the nota- 
tion ¥^:^{A) = f^FP{A)ep{dx)). 

Theorem 12.11 in combination with Theorem 17.5 of Sato |3QJ leads to the following result 
about the infinite time absolute ruin probability. We refer to Sato |30) for background on 
self-decomposable random variables. 

Corollary 2.2. 

(i) // log{y)i'{dy) < 00, then, under Fq ^\ Xt converges in distribution, as t ^ 00, to a 
positive self-decomposable random variable X^o and 

rP (to < 00) = Pj-''' (Xoo > X) . 

Moreover, for any f3 > 0, we have 

p5(ro = oo) = exp(-^'Md.). 

(ii) // /-|°° log{y)v{dy) = 00, then, for any x > 0, 

PW (to < 00) = 1. 

It is interesting to note that for the risk process without interest (r = 0), ruin is certain when 
the safety loading is negative that is whenever c < j-^ yz/(dy), whereas for r > the premium 
rate does not have any infiuence on whether ruin is certain or not. 
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Before stating the next two theorems we need to introduce a httle further notation. Let 
^ = iog{y)u{dy)=oo} ^^'^ define the measure W on [0, oo) via its Laplace transform as follows 

(2.2) / e-^^VF(dx) = e"^« " /3 > ^ 

Jo 

Noting, by an integration by parts, that 

ru Jq rx 

it follows readily that 

W{dx) = e«^/x(drE), 

where /i is the law of a positive, self-decomposable random variable with infinite Levy mea- 
sure '^^^'""•' dx (cf. [30, Corollary 15.11]). It follows that the function W{x) := W[0,x] 

is well-defined, increasing and by [30^ Theorem 27.4] that W{x) is continuous. Note that if 
log(y)z^(dy) < oo, then W is simply the distribution function of the random variable X^o 
from CoroUarv 12.21 We extend W to the whole real line by setting W{x) = on (— oo,0). Fur- 
ther, for all n < u{0,oo)/r, the n-th derivative of W, denoted by W^'^\ exists and is continuous 
on (-00,00) (cf. [30, Theorem 28.4]). 

The next theorem concerns a discrete spectral type representation of the transition distribu- 
te— r) 

tion of the process X under Pq . We stress that the spectral theory for self-adjoint operators 
in an Hilbert space structure is well established. In particular, McKean [23] discusses in details 
the nature of the spectrum of the semigroup of linear diffusions. However, this general theory 
does not apply here since we are dealing with non self-adjoint operators. In this context, there 
are little examples in the literature where such a spectral decomposition has been given. One 
notable exception is the paper of Ogura [26] • Therein, the author provides conditions under 
which the semigroup of continuous state branching processes with immigration (for short CBI) 
admits a discrete or continuous spectral representation. Furthermore, in a very elegant fashion, 
he manages to characterize through Laplace transform, the eigenmeasure and eigenfunctions 
associated to the semigroup. To be more precise, in [261 Theorem 3.1], Ogura shows, under 
some conditions, that for some tQ > the semigroup Pt{x,dy) of a CBI process satisfies the 
discrete spectral representation 

00 

(2.3) Pt{x, dy) = r/„(x)C„(dy)e-^"*, t>to,x,y> 0, 

n=0 

where rjn, resp. Cn, are eigenfunctions, resp. eigenmeasures, of Pt corresponding to the eigenvalue 
g-A„t ^[ii^ \^ > Q Although the process X under Pi "^^ does belong to the class of CBI processes, 
they were excluded in [26]. Moreover, the methodology of |26j does not extend to our case; in 
particular the suggestion made in footnote 3) on p. 309 of [26] does not lead to the right direction. 
In Theorem 12 . 31 below we are able to give, under some technical conditions on the Levy measure, 
a spectral representation similar to (]2.3p for the semigroup of this process when x = 0. We refer 
to Chapter 2.7 of Bingham et al. [5J for the definition of a quasi-monotone function and remark 
that the theorem remains valid if quasi-monotone is replaced by ultimately monotone. Recall 
that a function / : (0, 00) — )• (0, 00) is ultimately monotone at infinity if there exists j4 > such 
that / is monotone on (A, 00) and that / is slowly varying at infinity if lim2,._j.oo — 1 all 
t > 0. 
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Theorem 2.3. Assume that the Levy measure v satisfies 
i) e-^^jv(dx) < oo for all /3 G R, 

a) L'{x,oo) = i{l/x)x^°' with < a < 1 and £ a quasi-monotone slowly varying function at 
infinity. 

Then for any t > t^ = — ^ log cos (^) and x G M, the transition distribution of the process X 
under is given by 

oo 
n=0 

where 



Consequently, for any x > 0, 

oo 

(2.4) pM(ro>t) = ^/.„Ty(")(x)e-™*, t > t„. 

n=0 



v=0 



Moreover, for any n, W^^^^\y)dy is an eigenmeasure of Fx ^\Xt G dy) corresponding to the 
eigenvalue e"''"* in the sense that 

%"'\Xt G dy)W^''+^\x)dx = e-™*T^("+i)(y)dy, y G M. 

Our final theorem concerns the two-sided exit problem for the OU type process. More specif- 
ically, we compute the Laplace transform of the stopping time 

T+ = inf{t > : Xt > a}, 

on the event that X exits the interval [0, a] at o. When X is a spectrally negative Levy process, 
this quantity is given by fq{x)/fq{a) with q > being the parameter of the Laplace transform 
and where fq (resp. /o) is the so-called g-scale function (resp. scale function) of the spectrally 
negative Levy process, cf. |20l Theorem 8.1]. In Theorem 12.41 we derive a similar expression for 
the case where X is the OU-type risk process and show that the 'corresponding g-scale function' 
is given in terms of the fractional integral of the function W, which we denote by Wq. Hence if 
one has an explicit expression for W, one automatically gets Wq for g > in closed-form; this 
is in contrast to the case of spectrally negative Levy processes, see the discussion on p. 1674 of 
Kyprianou and Rivero [22] . 



Theorem 2.4. For q > 0, define the function Wq : {—oo, oo) — [0, oo) by Wq{x) = on (— oo, 0) 
and on [0, oo) by the Riemann-Liouville fractional integral of order q ofW, i.e. 

(2.5) Wq{x) = ^^ j\x-yY-^W{y)dy, x>Q. 

Then, for any q > 0, x < a and a > 0, we have 

^ {Ta <TOi ^ 



6 -^{r+<oo} 



Wq/r{x) 
Wq/r{a) 



We point out that Theorem 12.41 simultaneously gives the first-passage time above a of the 
OU-type risk process. For a process X defined by ()1.2p with Z a general Levy process with 
no negative jumps and with r < 0, Hadjiev [18] (under an extra assumption) and Novikov [25j 
provided an explicit expression for the Laplace transform of the first-passage time of X above a 
fixed level. Although not considered in these papers, it can be checked that their methodology 
and expression extend to the case where r > and Z a Levy process with no downward jumps 
and paths of unbounded variation. However, when r > and Z is a subordinator (that is, the 
OU-type risk process), the situation is completely different due to the possibility of absolute 
ruin. We remark here that we cannot expect to find a more explicit expression like the one 
given in and [25] for the r < case, since this would lead to an explicit expression for the 
distribution function of any positive self-decomposable random variable. 



3. Examples 

Before studying some examples let us recall the definition of a transformation recently in- 
troduced by Kyprianou and Patie [21j . which will be helpful in our context. For any Laplace 
exponent of a subordinator, that is of the form (jl.ip . we write for any 7 > 

Then, in [21], it is proved that Ty(/) is the Laplace exponent of the driftless subordinator with 
Levy measure i^-y(dx) = e~'^^(z^(dx) -|- 7z^(x, oo)dx). In particular, z^^(x,oo) = e~^^i'{x,oo). 
Obviously, 7o0 = We assume throughout this section that log(y)z^(dy) < 00 and therefore 
the quantity (friP) defined by 



is well-defined for /? > 0. It is easily seen, after a change of variable, that 

^{7)(^) ,= i f l^ds = Mf^ + 7)- Ml)- 

r Jo s 

In other words, the action of the mapping 7^ on the backward Laplace exponent is equivalent 
to the action of the Esscher transform on the Laplace exponent of the limiting distribution. In 
particular, if W{x; 7) stands for the limiting distribution function associated to the backward 
Laplace exponent Ty(j), then we have the following simple relationship 

(3.1) W'{x; 7) = e-'^^^^^-^'^W'ix; 0), x > 0. 

(Note that by [301 Theorem 28.4], W is always differentiable on (0,oo).) 

Next we detail some examples where the function W and/or the distribution of the absolute 
ruin time can be given in closed form. For more examples of cases where W is explicit, we refer 
to Iksanov and Jurek [19]. We end this section by giving an example that illustrates Theorem 
[231 

3.1. The compound Poisson case with exponential jumps. Let us start with the case 
where for x > 0, 

iy{x,oo) = r]e~^^ , 6,7] > 0. 

6 



Hence the subordinator Z is a compound Poisson process with exponentially distributed jumps. 
We have 

V..(/3) = ^log(l + /3/5) 
r 

and so the Laplace transform of X under Fi is easily computed via (j4.2p . By inverting this 

^(—r) 

Laplace transform, Perry et al. |29] shows that the transition function of X under Pq admits 
the following form 

^'^-'-\Xt G dy) = e-^'6o{dy) + ^e-''*(e'^* - l)e-^^Fi (l - ^; 2; 5{1 - e'')y) dy, y>0 

r \ r J 

where 5q stands for the dirac point mass at and iFi is the confluent hyper geometric function, 

see e.g. Lebedev |23j, Section 9.9]. Hence, from Theorem 12. H we deduce that 



M(ro>t) = e-''* 



1 + ^(e'^* - 1) £ e-^^Fi (l - ^; 2; 5(1 - e^')y) dy) 



From Corollary 12.21 it follows that ¥x\tq = oo) (or equivalently W{x)) equals the gamma 
distribution with shape parameter rj/r and scale parameter 1/5; this fact was first established 
by Gerber [T5] . 



3.2. The compound Poisson process with Linnik distribution. We consider an example 
found in Iksanov and Jurek |19) . Assume that the Levy measure of Z is given by 



v{x, oo) = riEa{—5x°'), x > 0, 

where ry, 5>0,0<a<l and 

oo 

fc=0 ^ ' 

is the Mittag-Leffler function and F is the gamma function. Hence the subordinator Z is a 
compound Poisson process with arrival rate tj and with jumps distributed according to a positive 
Linnik distribution. We have = -^fp^ and thus 

(^,(/3) = ^log(l + /375). 
ra 

We point out that in the case a = 1, the Linnik distribution boils down to the exponential 
distribution and hence this example can be seen as a generalization of the previous. Next, we 
deduce, from the identity (j4.2p below, that 

V 



ar 



Denote for k > by ^W{x) the increasing function on [0, oo) characterized through its Laplace 
transform which is given by 



e 



-I3x 



^w{dx) = (rV" + 1) " , /3 > 0. 



Note that jlW{x) = rP{To = oo). By means of the binomial formula, we get, for any /3° > 1, 

ar 

(1 + rV")-'^ = f;(-i)"£i^±^(rV")-("+'^\ 

^-^ nil [K) 



n=0 
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and hence we obtain by a term-by-term inversion 

\K)T{a{n + K))n! 



71=0 



""^ I (a, ok) 



where i^i is a Wright hypergeometric function, see e.g. Braaksma [TJ Chap. 12]. Note that for 
a = 1, 

recovering the exponential jumps case studied above. Next, observing that for /3 big enough by 
the binomial theorem, 



(5~i/3°e-"'^* + 1 V"" f-^ l-e""*^* ~ ~ 



ar 



5-1/3" + 1 J V <5-ie-°^*^" + 1, 



=e-''* 



^ T(^)nl (1 + J-ie-°^«/3")"' 
we get, by means of Laplace transform inversion and Theorem 12. H 

(3.2) pW(ro >t)= e~^' ( 1 + e^* ^^^^^^(1 - e— *)\T^(e^*x) ] . 

\ n=l Var/ / 



Note that, for any k > 0, the power series 

oo 

E 



n! 

n.=0 



is analytic in the disc of radius 1. Since for any i > 0, |1 — e"^*! < 1 and, clearly nW{y) < 1 for 
any y > and n = 1, 2 . . ., we deduce that the series on the right-hand side of (j3.2p is uniformly 
convergent on x R+. Note also that the identity p.2p simplifies when — is an integer. In 
particular, when r] = ar, the finite-time ruin probability is simply a weighted sum of 1 and the 
infinite-time ruin probability, i.e. 

pM(ro >t)= e-*?* + (1 - e-''*)iVF(x). 

Moreover, when a = 1, we have 

°° r ( j_ R\ ^ V ( -L KIT- 

^1 r(2)n! Jo r(2)n! T{n) 

=S(e^t _ 1) £ e-Se^'y^F, (l + ^; 2; (e^* - l)6y) dy 

=5(e^* - 1) £ e-'y^Fi (l - ^; 2; (1 - e^')6y^ dy 

where the last equality follows from the identity e-^^iFi (a; b; y) = iFi {b — a; b; —y), see, after 
correcting the obvious typo, formula 9.11.2 in [23j. This is consistent with the previous example. 



We finally note that when < a < 1, we can use (j3.ip to get a new example where the infinite- 
time ruin probability is explicit. 



3.3. The stable subordinator. We assume that Z under P is a stable subordinator of index 
< Q < 1. Its Laplace exponent takes the form </>(/3) = cos(7i-a/2) ' — 0- It is easy to verify that 
in this case ifriP) = cos(-Ka/2)ar ' — 0- O'^^ check by taking Laplace transforms and using 
the scaling property of Z that, 

%^-''\Xt G dy) = P {v^/'^{t)Zi E dy] where vit) = llLlJ!!. 

V / ar 

Hence, by writing Pa for the distribution of a positive stable random variable of index a, we 

deduce in combination with Theorem 12.11 that 

P«(ro>t) = P«(x^;-V°(t)) 

and by taking the limit as t — )• cxd, we get from Corollarv 12.21 

pM(ro = oo) = P,((ar)i/-x). 

We mention that a series representation of stable densities can be found in Sato [31]. Additional 
interesting representations has also been derived by Schneider [35j. In particular, in the case 
a = 1/2, we recall that the stable distribution admits the well known simple expression 

(3.3) Py,{x)= r ^y-^l\-y^^y)Ay. 

Jo v27r 

We consider now the image of (p by the mapping 7^ which gives 

cos(7ra/2j 

This Laplace exponent is a specific instance of the family of characteristic exponents of truncated 
Levy processes constructed by Boyarchenko and Levendorskii [6], also called tempered stable 
processes. From the identity (13. ip . we deduce that 





3.4. Example from Theorem 12.31 We now give an example where we use the representa- 
tion (|2.4p for the finite-time absolute survival probability in Theorem 12.31 Assume z^(x,oo) = 
£(l/x)x~" with < a < 1 and 

£{x) = CI|^>^}, where C,A>0. 

Then is a Levy measure that satisfies all the conditions of Theorem 12.31 Note that this Levy 
measure corresponds to a subordinator with stable-like jumps near the origin and no jumps 
larger than A. This model with claims bounded in size, might for instance be used in the case 
where the insurance company has excess-of-loss reinsurance. A straightforward computation 
involving integration by parts, shows that 

Vr{(3) = ^ (rr(l - a;0,/3A) - A"" (l - e"''^)) , 

where T{a; zq, zi) = J^^ e~^t'^~^dt is the incomplete gamma function. We want to check how 
well the spectral representation (j2.4p performs. Since we do not know an explicit expression 
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for W, we proceed by numerical inversion. Let N G For n = 0,l,2...,A^ — 1 we compute 
^(n+i) inverting the right-hand side of (j4.9p below using the (inverse) fast Fourier transform. 
We then numerically integrate W' to get W and consequently we get an approximation for the 
right-hand side of (j2.4p where the upper boundary oo is replaced by A^. We remark here that 
as n increases, it gets more difficult to obtain accurate approximations for T^^""*"^) since these 
functions oscillate more and more as n increases. In particular, for n > 0, has n distinct 

zeros, cf. Theorem 5.1]. In order to get an approximation for the left-hand side of (12. 4p . 
note that from (14. 2p , it follows that the Fourier transform of the transition density of Xt under 



ii-r) 



is given by 



(3.4) 



E, 



(-r) 



u G 



t > 0, 



and by numerical Fourier inversion and integration, one gets (cf. Theorem l2.ip an approximation 
of the finite-time absolute survival probability at a fixed time t. Hence we obtain numerical 
approximations of the truncation error corresponding to (12. 4p which we define by 



max 
=0,1,. ..,M 



N 



hi 



To 



>t)- J]/z„VF(")(/i^ 



^ e 



-rnt 



n=0 



For our example we take r = 0.2, a = 0.5, C = A = 1, h = 0.2 and M = 125. All the calculations 
are done in Matlab. In Tabled! numerical approximations of the truncation error are displayed 
for different values of and t. One sees that convergence takes place for t = 7, 10, 15 as grows 
and that for a fixed A^, the error becomes smaller as t grows. We point out that the numbers 
in the table, besides the truncation error, also consists of integration error due to the Fourier 
inversion. In particular, the integration error will dominate the truncation error for high values 
of A^ (and low values of t) due to the highly oscillating nature of W^"^ for large values of n; this 
might explain the large value in the table for A^ = 16 and t = 5. The spectral representation 
performs badly for t = 3, but note that 3 < « 3.46 and hence this case is not covered by 



jn=0 



„VF(")(a;)e-™* to {tq > t) for t = 7 is graphically 



Theorem 12.31 The convergence of 
displayed in Figure 1. 

We remark that the benefit of computing finite-time absolute survival/ruin probabilities via 
(j2.4p is that one can quickly calculate these probabilities for a whole range of time points t, 
whereas if ones uses the method via (j3.4p . one has to perform a separate Fourier inversion for 
each t. Moreover, by using the spectral representation, one can take advantage of (j3.ip to quickly 
obtain the distribution of the absolute ruin time in the case where the tail of the Levy measure 
equals e~'^^i'{x, oo) with 7 > 0. 
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N 


3 


5 


7 


10 


15 





0.905 


0.718 


0.526 


0.312 


0.130 


1 


0.768 


0.461 


0.244 


0.091 


0.027 


2 


1.283 


0.453 


0.139 


0.025 


0.020 


3 


1.587 


0.385 


0.090 


0.024 


0.021 


4 


2.424 


0.426 


0.088 


0.025 


0.021 


6 


4.080 


0.349 


0.039 


0.022 


0.021 


9 


9.320 


0.237 


0.033 


0.022 


0.021 


12 


22.508 


0.167 


0.031 


0.022 


0.021 


16 


582.088 


0.887 


0.030 


0.022 


0.021 



Table 1. Numerical approximations of ejsi^t for different values of N and t. 






N=1 




N=3 








N=6 

N=oo 


y If 


■ N=0 



-0.2 L 



Figure 1. Graph of x ^ YZ=q ^nW^'^^x)^-'''^ for t = 7 and N = 0,l,3,6,oo. 



4. Proofs 

Proof of Theorem \2.1\ First, by means of a change of variables and the duality lemma for 
Levy processes, see e.g. [H p. 45], we obtain 



-r{t~u) 



d{Zt^u - Zt) G dy 



(4.1) 



^*e-'^(*-")dZ„Gdy 
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Let X > 0. Now using that the ruin time is absolute, ()1.2p and ()4.ip . 



e-^MZ, < -X 



^-'■^ {X, > x) , 



which gives the first assertion. The second follows directly from the well-known expression for 
the Laplace transform of Xt under Pi which reads, 



(4.2) 



exp ( -/3xe-'^* + J (/)(Ae-'^")ds ) , /3, t > 0, x G M, 



see e.g. Sato 



□ 



Denote by p/*^^ (resp. the semigroup of the process X under P^^^) (resp. p(-'')), i.e. for 

any positive measurable function / on R, 

P^^'-^fix)=KP[fiXt)], xGR,t>0. 

For the proof of Theorem [23] we need the following lemma; note that part (ii) is a (weak) duality 
result between the two semigroups. 

Lemma 4.1. 

(i) For any t > and x G M, 

(4.3) iXt,¥l) = ie^'x + Xt,ro), 

where = stands for the identity in distribution. 

(ii) For any two positive measurable functions f and g, 

(4.4) = e-^*(/,P^'-)5), 

where {f,g) = f{x)g{x)dx. 

(iii) For all a>0, (Xt € [-a, a] for all t>0) = 0. 

12 



Proof. Item (i) follows directly from (jl.2p . Next, using successively Tonelli, ()4.3p . a change of 
variables, (j4.ip and (|4.3p again, we get 

{Pt^f,9) = [ EP[fiXtMx)dx = E'^^^\ [ fiXt + e^'x)gix)dx 



(r) 



-rt 



-rt 



-rt 



-^'fiz)g{iz-Xt)e-^')dz 
t 





/(^)IE 
/(^)IE 



9 ze 



-rt 



9 e 



-rt 



z+ j\'-'d{-zs)y^ 



dz 



/(z)Ei-^)b(Xi)]dz 



=e-^\f,Pt^9), 

yielding (ii). Finally, choosing f{x) = l^x£[-a,a\} a-^d g{x) = l{x>o} with a > in (j4.4p leads to 

j^'S'^^\Xt e [-a,a])dx = 2ae~''* and now by applying standard argument from measure theory 
(cf. Exercise 3.1.12 of [3B]) and Fatou's lemma, we get 

/■OO f'OO 

/ p(r)(liminf{Xt G [-a,a]})da; < / liminf p(r)(Xt G [-a, a])dx 
Jo Jo 



< lim inf 



P{Xt G [-a,a])dx = 0. 



Hence F^x\Xt £ [—a, a] for alH > 0) = and the proof is completed. 



□ 



Proof of Theorem 2.4\ The first equality follows directly from Lemma I4.1l fiii) and (j2.ip . In 

order to prove the second equality, we first show that for all g > 0, 



(4.5) 



Let P > 0, q > and recall that ^ = ^{f^ iog(y)u{dy)=oo}- From the identity (j2.2p and observing 
the convolution in (|2.5p . we get 



(4.6) 



-'^''Wg{x)dx = l3-i'h-^i ^'^^ 



On the other hand, using successively (j4.4p . (j4.2p and (j4.6p . we obtain, denoting e/3(x) = e 
a; G M, that the Laplace transform of the left-hand side of ()4.5p is given by 

e-'?*(e^,p(^V,/,) = e-('?+^)*(P^'-)e^,Ty,/,) 



Hence by uniqueness of the Laplace transform and the continuity of Wq, ()4.5p follows. 
Now with the aid of (j4.5p we can use the Dynkin formula [11, Theorem 12.4], to derive 



(4.7) 
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Since the mapping Wg is increasing and continuous, the proof is completed by an argument of 

dominated convergence and the property that under = a on {t^ < 00} which follows 

by absence of positive jumps. □ 

Proof of Theorem \2.3[ Denote ipr{v) = Jq ^^du, v £ C. By assumption (i), ipr and conse- 
quently exp{{pr{-)) are entire functions (see e.g. [30l Section 25]) and hence we can write 

(4.8) e^^(^) = ^^^t;", veC. 

n=0 

By assumption (ii), z^(0,oo) = 00 and thus from Sato and Yamazato [32', Equation (2.20)], we 
get 

1 f °° 

(4.9) [y) = — (in)"e-^-(*")+'"f du, n = 0, 1, 2, . . . . 
From the above identity, we deduce that 

(4.10) MnH^("+')(y)t;" = Mn(ini')"e-'^^(^")+'"^d« 

and now we want to use Fubini to switch the sum and the integral. We have for w > by (j4.8p . 
J]] L„(iuf)"e-^'-(i")+^"S' du<V2 / ^A^„(|u|z;)"e^^(-^'-('"))du 

poo 

=V2 / e^'-(l"l^)e-^'(^'-('"))d7/. 
J —00 

Using Fubini, we have 

rlne((^r(i'u)) = / (1 — cos(ux)) dx 



/ / xsin(xs)ds dx 

Jo Jo X 

s\n{xs)i'{x^ oo)dxds. 



JO 

u roo 



10 JO 

Because of condition (ii), it follows by ^ Equation (4.3.8)] that 

j sin(xs)i^(x, oo)dx ~ s"~"'^£(s)r(l — a) cos ^— ^ , as s — )• 00 

and hence by Karamata's Theorem [SJ Theorem 1.5.11], 

^ , ,. .r(l — a) /vrax 

rv\t{(pr{iu)) ~ u £(u) cos — , as u — > 00. 

a \ 2 / 



Here f(x) ~ g(x) as a: — ?• co stands for lim^j-^oo = 1- Similarly for > by Fubini, 

g{x) 

Jo X 

xe ds— ^- — -dx 



/ / e~^^i'{x,oo)dxds 
Jo Jo 



10 Jo 

and thus by Karamata's theorem and Karamata's Tauberian theorem [Sj, Theorem 1.7.1], 

riprVu) ^ t(u) u , as u — >• cxD. 

a 

In particular ripr{u) is regularly varying at infinity with parameter a and thus 

rifrivu) rifrivu) rifr{u) ^ 

lim pr- r = hm — — r = V . 

It follows that r[^t{ipriiu)) - ^r{vu)] ~ [cos (^) as u ^ oo if w° / cos (^). 

Assume now that cos (^) > and let 5 = ^. Then d\t{ipr{}u)) — ipr{vu) is regularly varying 
at infinity with parameter a and it follows by e.g. [5l Theorem 1.5.6(iii)] that there exists some 
^7 > 1 and A > such that for all u > U, 9^e((^,.(i?/)) — ipr{vu) > Au""'^ . Since yic{(pr{iu)) 
and ipr{vu) are even and continuous functions, we can now conclude that the right-hand side of 
(liTTTl is finite. 

This allows us to use Fubini in (|4.10p for < cos (^) to get with the aid of (j4.8p . 



n=0 ^^•^-O0n=0 



it/w)"e-^'-(^")+^"J'dtx 



1 

"2^ 



oo 

(pr(iTO)-'Pr(m)+ms/^^_ 



e 

oo 



The right-hand side of the previous identity is the inverse Fourier transform of the Fourier 
transform of '"^(O, dy) with v = e"''*. From earlier considerations, we see that this Fourier 
transform is integrable when t > 0. This implies in particular that for t > 0, '^^(0, dy) has a 

continuous density pl ^''(0, y), see e.g. the first line of Section 28 in Sato [30]. Since jyi y) is 
trivially integrable as well, we can now use the Fourier inversion theorem (see e.g. [141 Theorem 
8.26]) to conclude 

oo 

li~'\0,y) = j;^„I^("+i)(y)e-™*, t > log cos 

^-^ ra \ 2 

n=0 



The first identity follows then from the identity ()4.3p . 

To prove the second identity, noting that the right-hand side of ()4.1ip does not depend on y, 
we see that for any / : R — t- M measurable and integrable (i.e. \ f{y)\dy < oo), we can use 
Fubini to deduce 



oo „ ^ 

f{y)li"'^ (0, y)dy = V ;u„e-™* / /(2/)T4^("+i) {y)dy, t> log cos 



/ira 
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Choosing /(y) = l{y<2^} and applying Theorem 12.11 shows (|2.4p . The last part follows by taking 
Laplace transforms on both sides, hereby making use of Pubini, (j4.2p and (j2.2p . Note that by 
p. 190/191 of [30], W^^~^^\x) goes to zero as x goes to ±oo and hence W^"'~^^^ is bounded on 
(—00,00) which allows one to use Fubini. □ 
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